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Abstract. For the symbiotic branching model introduced in (EF04), it is shown that aging and intermittency 
exhibit different behaviour for negative, zero, and positive correlations. Our approach also provides an alternative, 
elementary proof and refinements of classical results concerning second moments of the parabolic Anderson model 
with Brownian potential. Some refinements to more general (also infinite range) kernels of recent aging results of 
(DD07) for interacting diffusions are given. 
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1. Introduction 

For the last three decades, equations of the type 

du{t,i)= ^ a(^, j)(u(t, j) - u{t, i))dt+ ^Kf{u{t, t)) dWt{i) (1.1) 

jei,'' 

have been studied intensively. Here, i G Z'', i > 0, k > 0, {a{i, are transition rates on Z'', and 

W = {Wt{i)\t > 0,i e Z''} is a familiy of independent Brownian motions. The following special cases with 
very different interpretations and different behaviour are quite common in the literature. 

Example 1. The (Wright-Fisher) stepping stone model from mathematical genetics: f{x) = a;(l — x). 

Example 2. The parabolic Anderson model (with Brownian potential) from mathematical physics: f{x) = x^. 

Example 3. The super random walk from pure probability theory: f{x) = x. 

Example 4. The critical (spatial) Ornstein-Uhlenbeck process: f{x) = 1. 

For the super random walk, k is the branching rate which in this case is time-space independent. In 
(DP98), a two type model based on two super random walks with time-space dependent branching was 
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introduced. The branching rate for one species is proportional to the value of the other species. More 
precisely, the authors considered 

du{t, i) = a{i, j){u{t, j) — u{t, i)) dt + a/ Ku{t, i)v{t, i) dWl{i), 
dv{t,i) = ^ a{i,j){v{t,j) ~ v{t,i))dt + Ku{t,i)v{t,i) dW^ (i), 

where now = {W^{i)\t >0,ie Z"^} and = {W^ii)\t > 0,i e Z'^} are families of independent Brow- 
nian motions. Solutions are called mutually catalytic branching processes. In the following years, properties 
of this model were well studied (see for instance (CKOO), (CDG04)). 

In this paper, we are interested in a variant of mutually catalytic branching, namely the symbiotic 
branching model introduced in (EF04) for continuous space. The same equations as for the mutually catalytic 
branching model are considered but additionally the driving noises are correlated in the following way: 

'■ i = j and n m, 
W"'{j))t = {t : ^ = j and n = m, (1.2) 
: otherwise, 

where G [— 1,1] is a correlation parameter. For g = and as well for general g there are basically two 
approaches to formalize the equations. In (DP98) under quite restrictive assumptions on the transition 
kernel (a(i, j))ijgz<i, existence of solutions was obtained in the space of tempered sequences. Since their 
assumptions in particular assume symmetry and exponential decay of («(«, ^gZ"* i already existence of 
solutions in cases we are interested in is not clear. This is why we stick to the setup of (CDG04), which as 
well is more popular for interacting diffusions. For the transition kernel (a(«, j))ijgz'' assume 

(Hi) 0<a(i,j)<oo, 

jezd 

(i/3) a{i,j) = a(0,z - j). 

Two main examples of interest are the following: 
Example 5. The discrete Laplacian is given by 



a{i,j) 



otherwise. 



Obviously, (Hi), {H2), {H3) arc fulfilled. Further, the one-dimensional Riemann walk (sec for instance 
(HugQS)) has transition rates 



with c normalising the total rate to 1. Here (i?i), (^^2), (^^3) are also fulfilled but in contrast to the discrete 
Laplacian the assumptions of (DP98) are not satisfied. 

To specify the state space we fix a positive, summable function a on Z'' satisfying 

^a(^)a(^,J) < A'a(j), Vj S 
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for some finite constant K. In (CDG04) a possible choice (see their Equation (4.13)) is given by 

oo 

igZd ,1 = 

where /3 is positive and summable, p^"'^ denote the n-step transition probabihties, and K > 1. This is needed 
to verify the generalized Mytnik self-duality which was introduced for the continuous space analogue model 
in Proposition 5 of (EF04). The duality for the discrete space is similar. For the duality in E in the special 
case g = see Lemma 4.1 of (CDG04). The state space is now defined by pairs of functions of the following 
Liggett-Spitzer space 

E={f:Z''^ M>o| ^ fUMj) < oo}. (1.3) 
The choice of a does not influence the results. 

Proposition 1.1. For uq,vo G E and g G [~1,1] there is a (weak) solution of the symbiotic branching 
model with almost surely continuous paths and state space E. For g E [—1,1) solutions are unique in law. 

The proof of Proposition 1.1 is standard. Existence can be proven by finite dimensional approximations 
as in (SS80). For g G (—1, 1) uniqueness follows from the generalized Mytnik self-duality as in (EF04). For 
g ~ ~1 uniqueness is true since moments increase slowly enough, and for g = 1 uniqueness of solutions is 
not known. 

For this work we restrict ourselves to homogeneous initial conditions 

uq ^vo = 1. 

This is not necessary but simplifies the notation a lot. 

The interesting feature of the symbiotic branching model is that it connects Examples 1-3 above. Being 
first established as a time-space inhomogeneous version of a pair of super random walks, the examples from 
above appear as special cases: g = obviously corresponds to the mutually catalytic branching model. The 
case g = —1 with the additional assumption uq + vq = 1 corresponds to the stepping stone model as can be 
seen as follows: Since in the perfectly negatively correlated case W^{i) = —W'^{i), the sum u + v solves a 
discrete heat equation and with the further assumption uq + vq = 1 stays constant for all time. Hence, for all 
t > 0, u{t, •) = 1 — v{t, ■) which shows that m is a solution of the stepping stone model with initial condition 
uq and w is a solution with initial condition vq. Finally, suppose w is a solution of the parabolic Anderson 
model, then, for g ^ 1, the pair {u,v) :— (w,w) is a solution of the symbiotic branching model with initial 
conditions mq = wq = Wo- 

The purpuse of this and the accompanying paper (BDE09) is to understand the nature of the symbiotic 
branching model better. How does the model depend on the correlation gl Are properties of the extremal 
cases g G {—1,0, 1} inherited by some regions of the parameters? Since the longtime behaviour of super 
random walk, stepping stone model, mutually catalytic branching model, and parabolic Anderson model are 
very different, one might guess that for varying g different regimes correspond to the different models. 

The focus of (BDE09) lies on the longtime behaviour in law (unifying the classical results for the stepping 
stone model, mutually catalytic branching model, and parabolic Anderson model) if (a(i, j))ijgz<i generates 
a recurrent Markov process, on (un)boundcdness of higher moments E[u(i, k)P] as i — > oo, and the wave 
speed for the continuous space analogue. It was shown that in the recurrent case E[u(i, k)P] is bounded in t 
if and only if 
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For the transient case the behaviour in g is open to a large extent. 

In contrast to (BDE09) the present paper focuses on second moments. Note that (1.4) implies that in 
the recurrent case second moments are bounded if and only if g < 0. This can also be seen and analysed in 
more detail using a moment-duality which will be explained in Section 3.1. Using this duality, we show how 
to reduce second moments of symbiotic branching processes to moment generating functions and Laplace 
transforms of local times of discrete space Markov processes, i.e. 

where k G M and Lt denotes the local time j^5f){Xs)ds in 0. For simple random walks the behaviour as 
t — > oo was partially analysed in (CM94) by analytic methods. Here, we present a simple new proof based 
on a renewal-type equation and Tauberian theorems. The simplicity of the proof has the advantage that no 
further assumptions on the Markov process (in particular no symmetry and no finite range assumptions) 
are needed. For any k > and k < the technique yields precise growth rates including all constants. 
As an application, intcrmittency and aging for symbiotic branching processes are established. 

The main results on intcrmittency and aging are collected in Section 2. In Section 3 we first establish 
representations for second moments of symbiotic branching processes (Section 3.1), then prove the results 
for exponential moments of local times (Section 3.2), and, finally, proofs of the main results (Section 3.3) 
are given. 

2. Results 

2.1. Intermittency 

The first property we address is intermittency (see for instance (CM94) for a discussion of the ideas). The 
p-th Lyapunov exponent is defined by 

j;{K,g)= lun ^\ogE[u{t,k)P] (2.1) 

if the limit exists (7p(K, g) analogously). Since in this work we only deal with second moments we further 
define 

^2''"{K,g)= lim -\og¥\u{t,k)v{t,k)]. 

In Lemma 3.2 we will see that 72'"(k, g) = 72 («^7 o) = 72 (^7 O) ^-iid hence we abbreviate 72. One says the 
system is intermittent (or weakly intermittent as recently in (FK09)) if 72 > 0. 

Intermittency for the parabolic Anderson model [g = 1) is a well-studied property (see (CM94) and 
(GdH07)). The existing proofs heavily depend on the linear structure of the system since they employ 
explicit solutions given as Feynman-Kac type representations. Such explicit solutions arc not known to exist 
for the symbiotic branching model. Hence, one might ask whether or not the results obtained for g ^ 1 can 
be transferred to some larger regime of correlation values. Indeed, this can be done. 

Let us first fix some notation. In the following [Xt ) denotes a continuous time Markov process with transition 
rates ia{i, and Ptihj) = ^[Xt = .jIXq = i]. Due to the moment-duality for symbiotic branching 

processes (see Lemma 3.8) the notation of symmetrization is needed. For two independent Markov processes 
(X^), (X^) with transition rates i))i_j(zz<i the symmetrization is defined as 

Xt=Xl~Xl (2.2) 

The transition rates of the symmetrization are given by 
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its transition probabilities arc denoted pt{i,j). Note that in the symmetric case pt(i,j) — P2t{i,j)- The 
Green function of {Xt) is denoted Gooihj) = Jq Pt{hj)dt and we abbreviate Goo — Goo(0,0). Further, 
we set Hooihj) = tpt{i,j)dt and abbreviate Hoo = ^?oo(0, 0). Analogously, we use Goo,Hoo for the 
symmetrization (Xt). 

Theorem 2.1 (Weak Intermittency for Symbiotic Branching). Let (ut,vt) be a solution of the symbiotic 
branching model with homogeneous initial conditions. Then {ut,Vt) is intermittent if and only if 

1 

Kg> 

In particular, there is no intermittency for non-positive g. 

The previous theorem suggests deviding symbiotic branching into two regimes in which the g > regime 
behaves like the parabolic Anderson model with respect to intermittency. 

Although after understanding the two-types particle moment-dual (Lemma 3.1) one sees that the problem 
can be treated as for gi = 1, we present a new proof. In (GdHO?) results of (CM94) for higher moments 
of the parabolic Anderson model were generalized to more general symmetric transitions than the discrete 
Laplacian. Here, in particular, complete results for the asymptotic behaviour (exponential and subcxponen- 
tial) of second moments of the parabolic Anderson model with arbitrary transitions (a(i, j))i.jgz<i are proven. 



In the course of the proofs we obtain the expression 



for the second Lyapunov exponent, where p~^ is the inverse of the Laplace transform of the return probabili- 
ties (sec Proposition 3.9). This expression, by Tauberian theorems, gives us the explicit asymptotic behaviour 
for the Lyapunov exponents as function of k. In the following, denotes strong asymptotic equivalence, i.e. 
hi ~ /i2 means lim/ii/ft,2 = 1- 

Proposition 2.2. Let {ut,vt) be a solution of the symbiotic branching model with homogeneous initial 
conditions. Then for ng > -tt— , the map 72 : [-7^ — , 00) K>o, k ^ 72 (k, g) has the following properties: 

i) 72 is strictly convex, 

ii) 72 (k) < i^g for all k, and -^^^ — > 1 for k — > 00, 
Hi) if pt{0, 0) ~ ct~" as t 00, a < 1, we have, as k s- 0, 



72 (t) 



f(cr(l -a)K£<)i/(i-") :0<a<l, 
\exp{-~{cKg)~^ + o{k~^)) : a = 1, 



iv) if pt (0,0) ^ ct " , as t ^ 00, a > 1, we have, as k \ > 0, 



72 (k) 



\^ cr(2-a) J 



1 < a< 2, 



^{ng - 1/Goo)(log l/{Kg - 1/Goo))-' : « = 2, 



G 

^^{hcg - l/Goo) : a > 2. 

Here, T denotes the Gamma function. 

Our approach has the further advantage that the growth rates in the critical and subcritical regimes 
follow directly: 

Proposition 2.3. Let {ut,vt) be a solution of the symbiotic branching model with homogeneous initial 
conditions. Ifpt{0,0) ~ ct~°', as t —>■ 00, then the following hold: 
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• Q > Q and a > 1 

i) If < then 



E[u{t,kf] ^— , ast^oo. 



a) If Kg = ^T— then, as t oo, 



cr(2°-Q)r(!i) ^" : 1 < a < 2, 
%^t : a > 2. 



-t^-" : a < 1, 
og{t) : a = 1, 
1 + kGoo : a > 1, 



E[u(t, kf] - <( Kclog(t) : a = 1, as t -> cx). 



g < 



:a<l, 

E[^(t, fc)^] ^ <! ^ 1 I . 1 •a>l 



e(i-e«Gao) 

For g ~ 1 (parabolic Anderson model), the subcxponential growth was partially analysed for finite range 
transitions in (DD07) (see their page 15). An example which was not included is for instance the Riemann 
walk defined in Example 5. Since in this case 

pt(0, 0) ~ ci"^/'^, as t ^ 00, (2.4) 

it serves as a convenient example for the above results which exhibits a precise recurrence/transience tran- 
sition at /3 = 1. Further, the simple random walk on Z'' is contained with 

pt(0,0) ~ c^"''/^ ast-^oo. (2.5) 

Combining the intcrmittency result with the extension of the results of (CKOO) given in (BDE09), we 
support the unstable behaviour of symbiotic branching for g > 0. It is quite standard (sec for instance 
(GM90)) that spatial processes being intermittent have a very local property: For large times the mass of 
the process is concentrated on few sites ("islands" ). Since for symbiotic branching, on each finite box solutions 
approach each constant configuration infinitely often, the islands do not stabilize. Since the diffusion function 
has the form Ku{t, k)v{t, k), we see that u will not produce high peaks if v is very small and vice versa. 
Hence, we suspect that u, v are concentrated on the same islands. Understanding the pathwise behaviour 
better is an ambitious task for the future. 



2.2. Aging 

Recently in (DD07) the concept of aging was discussed for certain classes of interacting diffusions. They say 
that aging (for linear test-functions) appears if the limit 

lim C0T[u{t,k),u{t + .s,k)] 

t,s — >oo 

depends on the choices of s and t. Aging does not appear if this is not the case. The main results of (DD07) 
were formulated with more general test-functions, though, restricted to finite range transitions. Differently, 
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the present technique is restricted to hnear test-functions but not to finite range transition. Our results 
suggest that neither finite range nor the hnearity of test-functions is crucial. Symmetry of the transitions is 
assumed as in (DD07). 

In (DD07) it is shown that no aging appears in the parabolic Anderson model (in our model g = 1) 
in any dimensions for the discrete Laplacian. Further, for the super random walk (in our model related to 
g = 0) it was shown that aging appears exactly in dimensions 1,2. This leads to the question if there are 
different phases for the symbiotic branching model. We show that the model exhibits three different regimes; 
an Anderson model like behaviour for g > 0, a super random walk like behaviour for g = 0, and a stepping 
stone model like behaviour for g < 0. The new case g < and Corollary 2.5 suggest that there are three 
regimes in which the most prominent examples fall. 

Theorem 2.4 (Aging for Symbiotic Branching). Let {ut,vt) be a solution of the symbiotic branching model 
with homogeneous initial conditions and a{i,j) — a{j,i). Then, ifpt(0,0) ~ ct~", as t ~> oo, the following 
is true. 

i) If g > 0, then no aging occurs for any a > 0. 
If g then 

- no aging occurs, for any a > 1, 

- lim(^s^oo,iog(s)/iog(t)=acor[w(t, fc),u(t-Hs,fc)] = (1 -a)+, for a = 1, 

- limt.s_oo.s=at cor[u(i, k), u{t -\- s, k)] = * ]-r| , for any a <l. 

(l+a)^ 

Hi) If g < 0, then 

- no aging occurs, for any a > 1, 

- limt,s^ooJog(s)/iog(t)=aCor[w(t, fc),M(i-Ks,fc)] = (1 -a)+, for a = 1, 

p(2r+a)"°(l-r)°"^dr 

- \imt.s^oo,s=at cor[u(i, k), u{t + s, k)] = " 2-°r(a)r(i-a) ' f'^'^ "-"^V " < 1- 

We emphasise that our proof of Theorem 2.4 can be applied to more general interacting diffusions. In 
particular, the examples from the introduction are included. For finite range transitions, i), ii), iii) of the 
following proposition were proven in (DD07) . 

Proposition 2.5. Consider solutions of (1.1) with homogeneous initial conditions. Then for 

i) f{u) ~ , aging appears as in Theorem 2.4 i), 
ii) < q;i < f{u) < ai, aging appears as in Theorem 2.4 H), 

iii) f{u) — u, aging appears as in Theorem 2.4 ii), 

iv) f{u) ~ m(1 — u), aging appears as in Theorem 2.4 Hi). 

In the cases in which aging occurs, the upper and lower limits are bounded by the stated values up to constants 
depending on f. 

3. Proofs 

3.1. Some Results on Symbiotic Branching 

We start with a discussion on how second moments of symbiotic branching processes can be reduced to 
exponential moments of local times. Let us first recall the two-types particle moment-dual introduced in 
Section 3.1 of (EF04). This will be used to calculate second moments explicitly. Since the dual Markov process 
is described formally in (EF04) we only sketch the pathwise behaviour. To find a suitable description of the 
mixed moment ¥l[u{t, k)'^v{t, k)^^], n + m particles are located at position k E 1/ . Each particle moves 
independently as a continuous time Markov process on Z"* with transition rates given by (a(i, i))ijgz<i. At 
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time 0, n particles have type 1, m particles have type 2. One particle of each pair changes its type when the 
time the two particles have spent at same sites with same type exceeds an independent exponential time 
with parameter k. Let 

= total collision time of all pairs of same type up to time 
hf ~ total collision time of all pairs of different type up to time t, 
ll{a) ~ number of particles of type 1 at site a at time t, 
l\{a) = number of particles of type 2 at site a at time t, 

Note that since there are only ti + many particles the infinite product is actually a finite product and 
hence well-defined. 

Lemma 3.1. Let {ut,Vt) be a solution of the symbiotic branching model with initial conditions uq,vo G E 
and Q £ [—1,1]. Then, for any k S Z"*, t >0, 

E[u(t,fc)"z;(i,/c)™] = E[{uo,voy'e'''-^^+B^t'^]. 

Though the proof for the moment-duality was given in (EF04) (see the proof of their Proposition 9) only 
for the discrete Laplacian we skip a proof. For general transitions j))i.j(=z'' the proof follows along the 
same lines. 

Note that for homogeneous initial conditions uq = vq = 1 the first factor in the expectation of the right-hand 
side equals 1. Lemma 3.1 in the special case g = 1, uq = vq = 1 was already stated in (CM94), reproven in 
(GdH07), and used to analyse the Lyapunov exponents of the parabolic Anderson model. 
For g 1, the difficulty of the dual process is based on the two stochastic effects: On the one hand, one has 
to deal with collision times of random walks which were analysed in (GdH07). Additionally, particles have 
types either 1 or 2 which change dynamically. 

Second moments are special since particles of different types do not change types anymore. Hence, when 
starting with two particles of same type there is precisely one event of changing types. This is used to obtain 
the following representation of second moments. 

Lemma 3.2. Let [ut, vt) be a solution of the symbiotic branching model with homogeneous initial conditions. 
Then, for any k G Z'', t > 0, 

E[u{t,k)v{t,k)] ^Eie'^^-^'], 

E[u{t,kr]^EHt,kr]^l' + f^''!\ = 

where Lt denotes the local time in of the symmetrization {Xt) defined in (2.2) started in 0. 



Proof. The first expression for the mixed second moment follows directly from Lemma 3.1: There are two 
particles which start with different types. Since pairs of particles of different types are never forced to change 
their types, they stay of different type for all time. Hence, Lf = 0, Lf = Lt for alH > and the assertion 
follows. 

For the second expression note that there is only one possible change of types. Starting with two particles 
of same types one of the types may change and the particles can not change their types again. Using 
independence of the particles and the exponential time we can make this explicit. Let Y be an exponential 
variable with parameter k, denote by X the law of the two independent Markov processes, and Lt their 
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X 


r f^' 


/ Ke~ 




L Jo 


KE[Lt] +E 


E 





collision local time. Integrating out the exponential variable leads to 

= E^^^[e«(^^+^^^?hy<iJ + E^^^[e'^(^^+^^^?)ly>Lj 
= E^ I Ke-''^e'"^+''^(-^'-^)dx] +E-^[e'*^^'E^[ly>Lj] 

J^* ne-'^e^ dx + E[e^^* e'''^*] : g ^ 0. 
This proves the assertion. □ 

Now we prepare for the proof of the aging result. 

Lemma 3.3. Let {ut, Vt) be a solution of the symbiotic branching model with homogeneous initial conditions 
and symmetric transitions {a{i, j))j_jizid ■ Then, for any k G Z'', i > 0, 

E[u{t, k)u{t + s, fc)] = 1 + K / P2r+s{k, k)E[e''^^'-''] dr 

Jo 

and similarly for v. 

Proof. The proof is only given for u since due to symmetry the same proof works for v. We first employ 
the standard pointwise representation of solutions 



{t,k) ^ 1+^ / pt-sii, k)^J Ku{s,i)v{s,i) dWl(i) 



(3.1) 



yielding 

E[u{t,k)u{t + s,k)] 



1 + E 



t+s 



Pt-r{i,k)^jKu{r,i)v{r,i)dW^{i) ^ / Pt+s-i(j, k)y^ ku{1, j)v{l, j) dWl{j) 



Further, since martingale increments are orthogonal this equals 



i+e[5: 

Now using independence of W^{i), W^{j) for i ^ j and Ito's isometry we continue the chain of equalities as 



Pt-r{i:k)^jKu{r,i)v{r,i)dW^{i) ^ / pt+s-i{j, k)^J nu{l,j)v{l,j) dWl{j) 



1+^E / pt^r{hk)pt+s-r{hk)Ku{r,i)v{r,i) 

%^Zd 



dr 



= 1 + / ^ Pt-r{i, k)pt+s^r{i, k)KE[u{r, i)v{r, i)] dr, 
^0 ^ezd 

where we were allowed to change the order of integration since all terms are non-negative. Using Lemma 
3.2, which in particular shows for homogeneous initial conditions that second moments do not depend on 
the spatial variable, symmetry of the transitions, and the Chapman-Kolmogorov equality, we finish with 

1+ / ^pt_.(fc,^W_.(^,A:)KE[e"^^-]rfr 

= 1 + K [ p2t+s-2r{k, fc)E[e"^^'-] dr = l + K [ P2r+s(.k, fc)E[e"'?^*-'-] dr. 
Jo Jo 

□ 
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Since wc are going to examine the second Lyapunov exponent 72 of solutions we give a simple argument 
which ensures existence of the exponent. 

Lemma 3.4. Let {ut, Vt) be a solution of the symbiotic branching model with homogeneous initial conditions. 
Then the Lyapunov exponent 72 exists. 

Proof. Note that to ensure existence of the limits 

lim - \ogE[u{t,kf] 

t—>oo t 

it suffices to show subadditivity of logE[u(t, k)'^]. Using Lemma 3.2 this is reduced to showing subadditivity 
of logE[e''^'^'], where Lf is the local time in of {Xt) started in 0. Thus, by conditioning on Xg, we get 

logE°[e''^-^'+=] = logE"[e'^^-^-'E^"=[e''«^'^*]] < logE"[e"^-^-'E°[e''«^'^*]] = logE°[e''^-^-i + logE°[e'"e'^*]. 

□ 



3.2. Exponential Moments of Local Times 



3.2. L Preliminaries 

In Lemma 3.2, we observed that in order to study second moments of symbiotic branching processes it 
suffices to study exponential moments of local times of the symmctrization (Xt)- We now take up this 
issue and discuss exponential moments of Lt in greater generality than needed for the symbiotic branching 
model. For the following let (Xt) be a time- homogeneous Markov process with countable state space S and 
transition kernel (a(i, j))i_j£s ■ In particular, the transition rates are not assumed to be symmetric. 
We start with a renewal-type equation for exponential moments of local times. 

Lemma 3.5. Let Lt be the local time of (Xt) in i G S for the process started in i. Then for k G ffi the 
following eguation holds: 



E[e''-^*] = 1 + K / prii,i)E[e''^'-'-]dr, t>0. 
Jo 

Proof. We use the exponential series to get 



(3.2) 



E[e 



KLt] 



E 



1 + E 



1 + E 



00 

n 



E 



,n— 
t 

(5, {XsJ . . . Si{XsJ dsn . . . dsi 



00 ft ft 
«=1 -^0 



(5i(XsJ . . . Si{Xs^Jdsn . . ■ ds2dsi 



The last step is justified by the fact that the function that is integrated is symmetric in all arguments and, 
thus, it suffices to integrate over a simplex. We can exchange sum and expectation and obtain that the last 
expression equals 

i + K '^""^ / • ■ ■ / ^[^^1 ^i,---,Xs„=i]dsn... ds2dsi. 

Jo J Si J Sn-l 



Due to the Markov property, the last expression equals 



ft ^ ft ft 

l + K Ps, {t, i) Y / ■ ■ • / ^[^s,-si = 0, . . . , Xs„^si = 0] ds„ 

•^0 Jsi J Sn-i 



. ds2dsi 
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and can be rewritten as 

l + ^y ps,{i,i) [y^n"-^ J ■J F[Xs^=0,...,Xs„=0]dsn...ds2jdsi. 

Using the same line of arguments backwards for the term in parenthesis, the assertion follows. □ 

Remark 3.6. A similar renewal-type equation as (3.2) can be shown with essentially the same proof for a 
discrete-time Markov process. It reads 

m 
n=0 

where Pn(i,i) is the return probability after n steps and L„ is the number of visits after n steps. 
Similar equations were obtained for symmetric Markov chains on Z'' in (MR92) using a completely different 
technique. Note that neither symmetry nor any structure of the set 5* is needed. The information on the 
geometry of S is completely encoded in pt(i, «). 

For the rest of this section we fix the Markov process (Xt), i G S*, and abbreviate 

f{t)^pt{i,i), ff(i)=E[e'=^']. 

The return probabilities pt{i,i) are always assumed to be strongly asymptotically equivalent to ct^", as 
t oo, for a > 0, as for instance for simple random walks on Z,'^ and the Riemann walk on Z. Further, 
/ is monotone, decreasing, positive with /(O) = 1, and g is monotone, increasing, positive with g(0) = 1. 
The Laplace transform for a function h on M>o is denoted by h and the convolution of two functions /, g is 
denoted f * g. In this notation Equation 3.2 reads 

= l + t>0. (3.3) 

Taking the Laplace transform of Equation (3.3) leads to 

5(A) = i + «;/(A)5(A), A>0. (3.4) 

Obviously, since / is bounded by 1. /(A) is always finite for all A > 0. A priori this is not true for g but if 
so, we obtain a useful representation from (3.4). 

Lemma 3.7. If g{\) < oo, then 

5(A) = . (3.5) 

A(l-«;/(A)) 

In the following we proceed in two steps. First, we use (3.3) to understand in which cases g{t) grows 
exponentially in t and discuss properties of the exponential growth rate. The following corresponence between 
exponential growth and finiteness of Laplace transforms holds (existence of the limit was proven in Lemma 
3.4): 

lim — logg{t) > c if and only if g{c) = oo. (3-6) 

t^oo t 

This observation is particularly important for the second step in which we discuss the behaviour of g{t) 
as t — > oo. In the cases in which g{t) grows subexponentially (3.6) implies that 5(A) < 00 for all A > 0. 
Hence, Lemma 3.7 can be used for all A > 0. The strategy in this case is the following: By assumption, 
the asymptotic behaviour of f(t) as t tends to infinity is known, namely ct~". Using Tauberian theorems 
the asymptotic behaviour of /(A) as A tends to zero can be deduced. By Lemma 3.7 this determines the 
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asymptotic behaviour of g{\) as A tends to zero. Using Tauberian theorems in the opposite direction, the 
asymptotic behaviom' of g{t) as t tends to infinity is obtained. 

To manage the transfer from the behaviour of / to / and back from g to g the following Tauberian 
theorems are used. They are taken from (BGT89) (see Theorem 1.7.6, Theorem 1.7.1, Corollary 8.1.7, and 
the considerations at the beginning of Section 8.1, §3). 

Lemma 3.8. Let h he a monotone function on M>o with h{0) = 1, then the following hold: 
i) If a < 1 and S E M, then h{t) ^ ct^"(logt)'' as t ^ oo if and only if 

as X ^ 0. 
ii) If h{t) ~ ct"^ as t ^ oo, then 

MA) -clog(l/A), 

as X ^ 0. 

Hi) If a > 1 and h{t) ^ ct~" as i — > oo, then I := h[t) dt < oo and 

rcT^^^., :l<a<2, 

I -h{X) ^ lcX\og{\) :a = 2, 

\xj^th{t)dt :a>2, 

as X ^ Q. 



3.2.2. Analysis of k> 0, Exponential Growth 

The main point of the analysis is the following representation of the exponential growth rate which follows 
directly from Lemma 3.5. 



Proposition 3.9. Let k > 0, then 



ProoL First, (3.6) implies that 



(k) lim ilogE[e«-^'] = /"^f-V 

t^oc t \ K/ 



(3.7) 



inf{A|g(A) < oo} = lim - logg{t). 

t^oo t 



Moreover, 



We are done if we can show 



/-(i)...f{A|/,A)<i}. 

{\\m < = {a|/(a) < 1}. 

First we show "C". Due to Lemma 3.5 we obtain g{X) = j + Kf{X)g{X) which imphes g{X) > Kg{X)f{X). 

Since g{X) < oo this shows that /(A) < i. 

Now we show "D". First, iterating (3.3) yields for fixed n 



git) = ^'(f" * 1) w + «"+^(/*("+^' * gm. 



i=0 
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Using f{t) <1 and g{t) = E[e''^*] < e''* yields 

K"+i(/*^"+iU.g)(t) = /" (s).g(i-s)ds 

< r fl^e'^C*--) ds < / e'^^*-^) ds ^ 0, 



as n — > C50. Hence, for fixed t > 

oo 
1=0 

Taking Laplace transforms we note that g{X) is finite if and only if the Laplace transform of the right-hand 
side is finite. However, using Fubini's theorem (note that only k > needs to be considered) we obtain 



(oo \ oo _j oo 

i=0 ) 4=0 i=0 



which is finite since we assumed k/(A) < 1. □ 

In particular, the previous result shows that understanding suffices to understand the exponential 
growth rates of £[6^^*]. This is not difficult due to the following observation: / is a strictly decreasing, convex 
function with /(O) = Goo- Hence, is a strictly decreasing, convex function with limA-,o /^^(A) = oo 
and /~^(A) = if and only if A > Goo- This implies that /^^(^) = precisely for A < ^r— . This and more 
properties of the exponential growth rate are collected in the following corollary. 

Corollary 3.10. Let k > and rin) = limt_»oo j logIE[e'*'^']. Then with Her '■= jj— the following hold: 

i) r{K) > and r{K) > if and only if k > Kcr, 
ii) the function k ^ r{K) is strictly convex for k > Kcr, 
Hi) r{K) < K for all k, and —^ I, as k — > oo, 
iv) if a <\, then Kcr = and, as k —> 0, 



r{K) 

v) if a> \, then Kcr > and, as k\ k. 



K—{cT{l~a))— :0<a<l, 
cxp(— (c/t)"^ + o(k~^)) : a = 1, 



(---c)^(%t^)°"^ :l<a<2, 

.(k-Hc)^ :a>2. 

Proof. Parts i) and ii) are proven as argued above the corollary. 
Since / < 1, the first part of iii) follows from 



-I /"OO /"OO 

- = f{r{K))= / e-''('')^/(x) < / e-''('')^dx 



n Jo Jo 

Continuity of / and /(O) = 1 imply that for e > there is xo{e) such that f{x) > 1 — e for x < xo{e). Hence, 

- = /(r(^))= / e--<^^^-fix)dx 



K 



> (1 - e) / e-'-^")^ dx^il- e)—-(l - e 



r{K)'- 
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Since r(K) oo for k — » (X) we obtain 



liniinflM > 1. 



The second part of iii) now follows since as well ^ < 1 for all k > 0. 

Finally, for iv) and v) note that the asymptotics of / for A arc known from Lemma 3.8. This translates 
to f-^{\) and hence to r{n) = f~^{^)- □ 

3.2.3. Analysis of k> 0, Suhexponential Growth 

So far, we have understood the behaviour of E[e'^^*] = g{t) as t — > c» for k > ^r— . In this case g{t) grows 
exponentially and the behaviour of the exponential rates in k could be analysed. We now come to the case 
K < ^T— . First, if Goo = oo, there is nothing to be done since the only appearing case is k = which yields 
g{t) = 1 for all t>{). Hence, we can stick to Goo < oo. 

Proposition 3.11. Let k > and k < Then, as t ^ oo, 

E[e'=-^'l 



1 - kGc 



Proof. Since Goo < oo we can apply part iii) of Lemma 3.8. Hence, /(A) Goo, as A ^ 0. As discussed 
above, since g{t) does not grow exponentially, g{X) > for all A > 0, and we can use Lemma 3.7. This 
implies 

5(A) ~ A-^ ^ 



(1-^^Goo)' 

as A 0. Going backwards with Lemma 3.8, part i), a = S = 0, the asymptotic of g follows. □ 
Proposition 3.12. Let k > and n ~ — !— . Then, as t oo, 

^ Kcr(2'-QV(") : 1 < " < 2, 

Proof. Since Goo < oo we can apply Lemma 3.8, part iii). Hence, /(A) ^ Goo, as A ^ 0. As discussed 
above, since g{t) does not grow exponentially, g{X) > for all A > 0, and we can use Lemma 3.7. Since 
kGoo = 1, the denominator (1 — k/(A)) appearing in Lemma 3.7 does not behave like a constant and we 
cannot apply part i) of Lemma 3.8 with a = 6 = 0. Instead we use Lemma 3.8, part iii), to obtain 

A^-" ; 1< a < 2, 



cr(2-a) 



A- Goo -./(A) ,^>2, 
as A ^ 0. This, by Lemma 3.8, part i), implies the assertion. □ 
3.2.4. Analysis of n < 

We now investigate Equation (3.2) for k < 0. 
Proposition 3.13. If k < 0, then, as t 00, 

1 1 



-Kcr(l-a)r(a) 

L J I log t —KC 

1 



< a < 1, 
a = 1, 
: a > 1. 
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Proof. First note that for k < 0, g{t) = E[e''^*] < 1 and hence for all A > 0, ^(A) < oo which validates the 
use of Lemma 3.7. This implies 



A(l - ^f(A)) /(A) 

as A — > 0. Using Lemma 3.8 in both directions returns the assertion. □ 

3.3. Proofs of the Main Results 

Proofs of Theorem 2.1 and Proposition 2.2. These follow from Lemma 3.2 and Corollary 3.10. □ 

Proof of Proposition 2.3. This follows from Lemma 3.2 and Propositions 3.11 and 3.12. □ 

Proof of Theorem 2.4. For this proof we always denote by ^ the strong asymptotics at infinity and we 
abbreviate pt = pt(fc, k). Lemma 3.3 implies 



cor[u(t, k), uit + s, k)] = 



Step 1, Q ~ 0; First, assume a > 1 which implies p2rdr < oo. Since 

ft+s 



Jo Js Js 

we obtain, independently of the choice of t and s, 

r /, I \ /, T \ i /n P2r+s dv s,i— *oo ^ 

cor[uft,fc),Mft + s,fc)] = , ^ 0. 

Y /o P2r dr /J^'* P2r dr 

Here, we used / « 5 if < liminf f / g < limsupj/g < cxd. We now come to the case a ~ 1, where we get 
^ P2r+sdr r^c {2r + s + l)-^dr = | log( ^^^^^ 

Therefore, we have 

coT[u{t, k),u{t + s, fc)] ^ ' (3.8) 

Vlog(2i+l)log(2(i + s) + l) 

For s — with a < 1 this expression behaves asymptotically as 

^£££=lL = l-a 
v/log(<) log(i) 

On the other hand, for s = t° with a > 1 the term in (3.8) behaves asymptotically as 

log(l) s,t-»oo , 

Vlog(2(t + t-) + 1) log(2(t + t-) + 1) 

Hence, for log(s)/ log(i) = a, we obtain cor[u(t. A;), u(t + s, /c)] ^ (1 — a)+. 
Now suppose a < 1. Then 



0. 



I p2r+sdrr^cf {2r + s + l)-"dr = 
Jo Jo 



c (2t + s + 1)1-° - js + iy 
2 1 - a 
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Therefore, we have 

(2i + s + 1)1-" - (s + 1)1 



cor[M(t, fc), uit + s, k)] 



V((2i + 1)1-" - 1) mt +s) + 1)1-" - 1) ■ 

For s ~ at this behaves asymptotically as 

(2 + a)i-" - fli-" _ (1 + a/2)i-" - (a/2)i-" 
V^2i-" (2(l + a))i-" " (l + a)(i-")/2 ' 

Step 2, Q < 0; Let us first consider a > 1. Since Ci < Ee^"^*"'' < C2 this case is exactly the same as = 0, 
a > 1. Now suppose a = 1. In this case we have by Proposition 3.13 

c f*"'^ 1 c Z"*-! 1 1 

P2r+s^e^'^^^~- dr / — -dr ^ / — ^ — —dr. 

^ -ngcJi (2r + 1 + s) log(t - r) -kqc 2(i - r) + 1 + s log(r) 

We use the scaling s = t°- with a < 1. Let < 6' < 1. The integral above can be split from e to Ot and 6t to 
t — 1. We treat the first integral and show that its order is less than (log(i))-i. First note that in the range 
of integration 

1 1 1 

< < 



2(t - e) + 1 + s ~ 2(t - r) + 1 + s " 2(t - St) + 1 + s ' 
Therefore, 



dr K, — I —dr 



2{t -r) + l + s log(r) t log(r) log(t) ' 

On the other hand, the second integral can be treated as follows. In its range of integration we have 
< ™ < ^rrkr. Therefore, 



log(t) ^ log(r) ^ log(et) 



1 1,1 r~' 1 , 11, f2{t-et) + s 

dr ~ - — -— / — dr = - — -— — log 



et 



2{t-r) + l + s log(r) log(t) Jg^ 2{t - r) + 1 + s log(i) 2 V 2 + s 



Thus, 



cor[.(t, k), uit + k)] ^ ^ ^ = 1 - a. 



Analogously, for case a > 1. Therefore, we get cor[u(t, k),u{t + s, k)] ~ (1 — a)+, whenever log(s)/ log(t) = a. 
For g < only a < 1 is left: Here, we have 



t 



(2r+l + s)-"(t-r)"-idr 



Jo -KQcT{a)T{l - a) Jq 

We set s — at. The integral can be rewritten as 

L 

(2ri + l + af)-"(t-ir)"-4(ir ~ / (2r + a)-"(l - r)"-idr. 



The same way one can see that 



*p2.E[e^''^'-'-]dr —4— /\2r)-"(l - r)"-idr. 
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Thus, 

(2r + a)-"(l -r)"-idr 



cor[u(t, /c), u{t + s, k)] 



/g (2r)-"(l -r)"-idr 



2-°'B{a,l-a) 2-"r(a)r(l - a) 

when s ~ at —^ oo. Here, B denotes the Beta function. 

Step 3, g > 0: The transient case (a > 1) with kq < ^r— has aheady appeared in the case g = for 
a > 1. E[e^''-^'] is bounded due to Proposition 3.11. Hence, no aging occurs. 

For Kg > 7-!—, and for a < 1, we proved in CoroUary 3.10 i) that Efe^'^-^'l grows exponentially. This implies 
that there is a A > depending on the growth rate, i.e. on g and k, (for example A = i-P^^(^) does the 
job) such that for all t > and < r < t we have 

Therefore, 

/ p2r+s^eB''^'—] dr < PsEie^""^'] f e~^'' dr < \ PsEie^""^']. 
Jo Jo ^ 

On the other hand, by the assumption pt ^ ct^°', the renewal-type equation of Lemma 3.5, and since 

e> 0, 

/ p2rE[e^'^-^*-'-] dr>c f p^E[e^«-^'-] dr = — (EieS^H - l) > c'Ele'^^H. 

Jo Jo en 

Putting these pieces together we obtain that 



c.nr[u{t k),u{t + .', k)]<r" ^\ ^ [ ^ r"p^ ' ^ ^ 



Note that, since g > 0, the term with the square root is bounded by 1. Since clearly ps tends to zero, the 
whole expression must tend to zero independently of how t, s — > 00. 
The only case left is a > 1 and Kg = -p^. First we consider 1 < a < 2. 



rE[e'^'^^*-'']dr ^ / (r + l)""(t - r)""! dr = / (r + l/t)-"(l - r)""! dr. 



P2r 



This expression tends to infinity for t 00. The rate is 



1/2 .1/2 

(r + l/t)-"(l - r)"-i dr w / (r + l/t)"" dr « (l/t)^-" = t"'^ . 





On the other hand. 



P2r 







+,E[e^'^^*-^-]dr ^ / (r + s)-"(t - r)""! = / (r + s/t)""(l - r)""! dr. 



This expression is bounded or tends to infinity (deponing on how t, s ^ 00). It is bounded by 

1/2 .1/2 

(r + s/t)-"(l -r)"-idr « / {r + s/ty dr < ci + C2{s/t)^-°' . 
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Putting these pieces together we obtain that 



The calculation is completely analogous in the cases a>2. □ 

Proof of Proposition 2.5. The pointwise representation of (3.1) is not restricted to the symbiotic branch- 
ing model but also holds for the interacting diffusions of this proposition. Hence, the same derivation as for 
the symbiotic branching model yields 



\l /„' vAk. k)m[Huit - r, k))\ dr /„'+' piJk, 4)E|/(»(f + - r, t))] * 

Part i) is contained in Theorem 2.1 since the parabolic Anderson model appears as special case g — 1. For 
part ii) we can estimate the expectations from above and below to obtain the same result as in Theorem 2.4 ii) 
except constants. The same is true for part iii) since the pointwise representation (3.1) implies E[w(t, k)] — 1. 
Finally, we could interprete part iv) as a submodel of the symbiotic branching model. Instead, we give a 
direct proof using use the coalescing particles dual of (Shi88). The dual process consists of two independent 
particles started in fc, performing transitions {a{i, in continuous time. After spending an exponential 

time Y with parameter k, independent of the particles, at same sites, the particles coalesce. We denote X 
the law of the particles and suppose uq = w € (0, 1). Then 

E[w(t k)'^] ^ ^ [tf^^^^^^^^^ non-coalcsccd particlcsj 

= E^><^[^«ly<iJ +E^><^[ii;2ly>iJ = w{l - E^[e-«^']) + iw^E^ie-''^'], 

where Lt denotes the collision time of the particles. Using 'E[u{t, k)] = w this yields E[/(u(t, k))] ~ K[u{t, k)] — 
E[u(t, fc)^] = {w — ij;^)E[e~''^*] and we can proceed as for the symbiotic branching model with g < 0. □ 
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